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Abstract. We study intermediate sums, interpolating between 
integrals and discrete sums, which were introduced by A. Barvi- 
nok [Computing the Ehrhart quasi-polynomial of a rational simplex, 
Math. Comp. 75 (2006), 1449-1466]. For a given semi-rational 
polytope p and a rational subspace L, we integrate a given polyno- 
mial function h over all lattice slices of the polytope p parallel to 
the subspace L and sum up the integrals. We first develop an al- 
gorithmic theory of parametric intermediate generating functions. 
Then we study the Ehrhart theory of these intermediate sums, 
that is, the dependence of the result as a function of a dilation of 
the polytope. We provide an algorithm to compute the resulting 
Ehrhart quasi-polynomials in the form of explicit step polynomi- 
als. These formulas are naturally valid for real (not just integer) 
dilations and thus provide a direct approach to real Ehrhart theory. 
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1. Introduction 

Let p be a rational polytope in \^ = M*^ and h{x) a polynomial 
function on V. A classical problem is to compute the sum of values of 
h{x) over the set of integral points of p, 

s{p,h)= J2 H^)- 

The sum S{p,h) has generalizations, the intermediate sums S^{p,h), 
where L OV is a rational vector subspace, introduced by Barvinok [6]. 
They interpolate between the discrete sum S{p,h) and the integral 
Jp h{x) dx. For a polytope p C and a polynomial h{x) 

^ Jpn{x+L) 

where the summation index x runs over the projected lattice in V/L. 
In other words, the polytope p is sliced along affine subspaces parallel 
to L through lattice points and the integrals of h over the slices are 
added up. For L = V, there is only one term and S'^(p, h) is just the 
integral of h{y) over p, while, for L = {0}, we recover the discrete sum 
S{p,h). As in the discrete powerful method is to consider the 

intermediate generating function 

^ Jfin(x+L) 

for 1^ e y*. If p is a polyhedron, not necessarily compact, the generating 
function still makes sense as a meromorphic function. By writing a 
polyhedron as the sum of its cones at vertices (Brion's theorem), we 
need only study the case where p is an affine cone. 

It is then natural to turn to the Ehrhart theory of the intermediate 
sums S^, that is, the study of the intermediate sum S^(tp, h) of a di- 
lated polytope tp as a function of the dilation parameter t. It turns out 
that, just like in the classical case, the Ehrhart function t ^^(tp, h) 
is a quasi-polynomial, that is, a function of the form 

d+M 

S\tp,h)=J2E^;^{p,h;t)t"^, (2) 

m=0 

where the coefficients E^{p, h;t) depend only on t mod q, where q is 
an integer such that gp has lattice vertices. 
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The main results of this article are: 

1. a polynomial time algorithm for the computation of the intermediate 
generating function of a simplicial affine cone, when the slicing space 
has fixed codimension, Theorem 24, 

2. a polynomial time algorithm for the computation of the weighted 
intermediate sum S^{p, h) of a simple polytope p (given by its ver- 
tices), and the corresponding Ehrhart quasipolynomial t S^{tp, h), 
when the slicing space has fixed codimension and the weight depends 
only on a fixed number of variables, or has fixed degree. Theorem 28. 

The main new feature of this article is the use of a decomposition of 
any simplicial cone into simplicial cones which have a face parallel to 
the given space L, the decomposition being given by a closed formula. 
Theorem 18. This decomposition is borrowed from [*)] where it plays 
a key role in the study of partition functions. We prove that this 
decomposition is done by a polynomial algorithm, if the codimension 
of L is fixed. 

Thus, the computation of the intermediate generating function of a 
simplicial affine cone is reduced to the case where L is parallel to a face 
of the cone, and thus the generating function factors into a discrete sum 
and an integral. This case has been already studied in [2]. The output 
of the algorithm is a "short formula" where both ^ and the vertex s 
of the cone appear as symbolic variables (see Examples 6, 7 and 26). 
Once a short formula for the generating function is available, we follow 
the method of [2] for the computation of the weighted Ehrhart quasi- 
polynomial of a simple polytope. The method uses the fact that the 
intermediate generating function decomposes into meromorphic terms 
of homogeneous ^-degree. The output is again a short formula where 
the dilating parameter t appears as a symbolic variable inside step 
polynomials, see Example 30. 

The results remains valid for non-negative real values of t. Thus 
our method provides a direct proof of a real Ehrhart theorem (Theo- 
rem 31). This extension of Ehrhart theory to real dilation factors has 
been studied by E. Linke [ ] for the classical case. It is even more 
natural for the intermediate valuations, as one of the cases is L = V. 
Then S^(tp, h) is the integral of h over the tp, the result being clearly 
a polynomial formula valid for any non-negative real value of t. 

The present manuscript together with the article [j] also provide the 
foundation for a future work, in which the weighted Barvinok patched 
sums \ {L)S^{p, h) for a simple polytope will be studied. 
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2. Notations and basic facts 

2.1. Rational and semi-rational polyhedra. We consider a ratio- 
nal vector space V of dimension d, that is to say a finite dimensional 
real vector space with a lattice denoted by A. We will need to con- 
sider subspaccs and quotient spaces of V, this is why wc cannot just 
let \^ = M"* and A = Z'^. An element v E V is called rational if G A 
for some non-zero integer n. The set of rational points in V is denoted 
by Vq. A subspace L of y is called rational if Ln A is a lattice in L. If 
L is a rational subspace, the image of A in y/L is a lattice in V/L, so 
that V/L is a rational vector space. The image of A in V/L is called 
the projected lattice. It is denoted by Ay/L- A rational space V, with 
lattice A, has a canonical Lebesgue measure dx — dmA{x), for which 
V/A has measure 1. 

A convex polyhedron p in K (we will simply say polyhedron) is, by def- 
inition, the intersection of a finite number of closed half spaces bounded 
by affine hyperplanes. If the hyperplanes are rational, we say that the 
polyhedron is rational. If the hyperplanes have rational directions, we 
say that the polyhedron is semi-rational. For instance, if p C is a 
rational polyhedron, t is a real number and s is any point in V, then 
the dilated polyhedron tp and the translated polyhedron s-\-p are semi- 
rational. Unless stated otherwise, a polyhedron will be assumed to be 
semi-rational. 

In this article, a cone is a polyhedral cone (with vertex 0) and an 
affine cone is a translated set s-|- c of a cone c. A cone c is called simpli- 
cial if it is generated by independent elements of V. A simplicial cone c 
is called unimodular if it is generated by independent integral vectors 
Vi, . . . , f fc such that {vi, . . . ,Vk\ can be completed to an integral basis 
of V . An affine cone a is called simplicial (resp. simplicial unimodular) 
if the associated cone is. 

A polytope p is a compact polyhedron. The set of vertices of p is 
denoted by V(p). For each vertex s, the cone of feasible directions at s 
is denoted by c^. 

We denote by J0,{V) the space of functions on V which is generated 
by indicator functions [q] of polyhedra q which contain lines. 

2.2. Valuations and generating functions. 

Definition 1. Let M be a vector space. An M -valued valuation $ is 
a map from the set of polyhedra p gV to the vector space M such that 
whenever the indicator functions [pi] of a family of polyhedra pi satisfy 
a linear relation ri [pi] — 0, then the elements $(pi) satisfy the same 
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relation 

i 

Definition 2. We denote by T-L{V*) the ring of holomorphic functions 
defined around G V* . We denote by Ai{V*) the ring of meromorphic 
functions defined around E V* and by AiiiV*) C A4(y*) the subring 
consisting of meromorphic functions 0(^) which can be written as a 
quotient of a holomorphic function by a product of linear forms. 

The intermediate generating functions of polyliedra wliich are stud- 
ied here are important examples of functions in Ai^CV*). 

Proposition 3. Let L C V be a rational subspace. There exists a 
unique valuation which to every semi-rational polyhedron p gV as- 
sociates a meromorphic function ^^(p) G M.{V*) so that the following 
properties hold: 

(a) If p contains a line, then ^^(p) = 0. 
(b) 

^^(P)(0= E / e^^'yUmUy), (3) 

for every ^ & V* such that the above sum converges. 
(c) For every point s G A + L, we have 

5^(. + p)(0=e«'^)5^(p)(e). 

When necessary we will indicate the lattice and use the notation 
5^(p,A). 

Remark 4. The fact that 5''^(p) is actually an element of A4i{V*) will 
be a consequence of the explicit computations of the next section. 

For L = {0}, we recover the valuation S given by 

s{pm = E ^^''^^ 

xepnA 

provided this sum is convergent. For L = V, 5'^(p)(^) is the integral 

I{pm= /e«'^Mx, 

if p is full dimensional, and 5'^(p) = otherwise. If p is compact, the 
meromorphic function S'^(p)(^) is actually regular at ^ = 0, and its 
value for ^ = is the valuation E]^±{p) considered by Barvinok [6]. 
The proof is entirely analogous to the cases L = {0} and L = V, see 
for instance Theorem 3.1 in [7] and we omit it. 
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Remark 5 (Rational vs. semi-rational polyhedra). Although [7] and 
other classical texts deal only with rational polyhedra, their proofs ex- 
tend easily to semi-rational ones. For instance, by summing geometric 
series, one obtains immediately the generating function of a unimodular 
cone with a real vertex. 

Here are some simple examples. For t G M we denote by {t} the 
fractional part of t, i.e., < {t} < 1 and t — {t} G Z, and by [t] the 
ceiling of t, i.e., \t] = t -\- {—t}. 

Example 6. For any t G M, we have 

^(t + M>o)(e) = ^=e*«^. 

More generahy, let c C M'' be a unimodular cone with primitive edge gener- 
ators vi, . . . , Vd- Let s be any point in M'^, with s = Yli=i ^i'^i- Then 



i=l ' ~ 1=1 ^ 



Example 7. Let c C be the first quadrant and L = Mei. We write 
s = {a,b) G M2. We have 



°° f+OO y CO ^ 

S'-is + c)(6, 6) = E y e^^"""^' dx = ^ e<^ j 



: !_ = e^'''^^- - 



1 - e?2 ^1 1 - e?2 ^1 

Example 8. Let p C be the triangle with vertices (0,0), (1,0), (0,1). 

— c 

6 



Let L = M(1,0). Straightforward calculations give S^{p) = ^t^, /(p)(0 



^^(i^ - ^17^)' and 5(p)(0 = 1 + e«i +e«2. These three functions are 
indeed analytic. 

A consequence of the valuation property is the following fundamental 
theorem. It follows from the Brion-Lawrence-Varchenko decomposi- 
tion of a polytope into the supporting cones at its vertices [8]. 

Theorem 9 (Brion's theorem). Let p be a polytope with set of vertices 
V(p). For each vertex s, let Cs be the cone of feasible directions at s. 
Let L ^ V be a rational subspace. Then 

S^ip)= J2 sHs + Cs). 

sev(p) 
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Figure 1. The projected lattice A{i} 

3. Intermediate generating function for a simplicial 

affine cone 

3.1. Where the shcing space is parallel to a face. Let c C V 

be a simplicial cone with integral generators {vj,j = l,...,d), and 
let s & V. Let us recall the expression for the intermediate generating 
function S^{s + c){^), when the slicing subspace L is generated by some 
of the fj's. 

For J C {1, . . . we denote by Lj the linear span of the vectors 
{vi,i G /). We denote the complement of I in {1,. . . ,d} by I'^. We 
have V = Lj (B Ljc For x & V we denote the components by 

X = Xj ~t~ Xjc. 

Thus we identify the quotient V/Lj with Ljc and we denote the pro- 
jected lattice by A/c c Ljc Note that Ljc fi A C A/c, but the inclusion 
is strict in general. 

Example 10. Let vi = (1, 0), V2 = (1, 2), I = {2}. The projected lattice A/c 
on Ljc = Rvi is Z^. See Figure L 

We denote by Cj the cone generated by the vectors Vj, for j E I and 
by bj the parallelepiped bj = ^jg/[0, l[fj. Similarly we denote by Cjc 
the cone generated by the vectors Vj, for j G P and bjc = If"^*- 
The projection of the cone c on V/Lj = Ljc identifies with tjc Note 
that a generator fj, i G I'^, may be non-primitive for the projected 
lattice A/c, even if it is primitive for A, as we see in the previous 
example. We write s = sj + sjc 

First let us recall from [ ] how the intermediate generating func- 
tion S^'{s + c, A) decomposes as a product. This already appeared in 
Example 7. 
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Proposition 11. [2]. The intermediate sum for the full cone s + c 
breaks up into the product 



+ c, A)(0 = S{si. + c,c,A,.)(0/(s, + c,,L, n A)(0. (4) 



We remark that in the above equation, the function S'(s/c+c/c, A/c)(^) 
is a meromorphic function on the space {Ljc)*. The integral I{si + 
C/jL/ n A)(^) is a meromorphic function on the space (Lj)*. We con- 
sider both as functions on V* through the decomposition V = Lj(BLjc. 

3.2. Cone decomposition with respect to a linear subspace. 

Let L be a linear subspace of V. Let c be a full dimensional simpli- 
cial cone in V with generators wi, . . . ,Wd- The main theorem of this 
section, Theorem 18, gives a signed decomposition of c into a family of 
full dimensional simplicial cones c^, where each has a face parallel 
to L. More precisely, this decomposition is a relation between indica- 
tor functions of cones, modulo the space CiV) of functions on V which 
is generated by indicator functions [q] of polyhedra q which contain 
lines. This decomposition is borrowed from [9]. The proof relies on 
the following simple lemma which has its own interest. The extreme 
cases, dimL = 1 or codimL = 1, of Theorem 18 follow easily from this 
lemma (Proposition 14). The general case is proven by induction on 
either dimL or codimL. 

Given vectors Ui,...,Up G M."^, we denote the cone Yl^=i^>o'^i by 
c(mi, . . . , Up) and its relative interior Yl^=i ^>oUi by c(mi, . . . , Up). 

Lemma 12. Let Ui, . . . ,Ud be a basis of V and let Ud+i = — J2k=i ^fc- 
For i = 1, . . . ,d + 1, define semi-open cones 

= c{Ui, + c(Mj+i, . . . , Ud+l). 

Then V is the disjoint union of the ti for i = 1, . . . ,d + 1. 

Proof. Let x = tiUi + . . . + tdUd- If t^ > for all k, then x G td+i- 
Otherwise, let i be the largest index with 1 < i < d, for which the 
minimum mmi<k<dtk is reached. Writing Ui = -J2i<k<d+i,kj^i'^k, we 



As tfc — tj > if 1 < j < (i and tj — ti > if i < j < d, and ti < 0, we 
see that x E ti. 

Let us show that ti fl tj is empty if z < j. Let x E ti, x = 
T.i<k<d+i,k^i^kUk. Writing Uj = -J2i<k<d+i,k^jUk, we obtain 



obtain 





l<k<d,k^i 



X = 





l<k<d,kf{i,j} 
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As tj > 0, it follows that x ^ ti. □ 

From the lemma, we deduce other interesting cone decompositions. 
The first one is the well-known stellar decomposition. 

Proposition 13 (Stellar decomposition). Let wi, . . . ,Wd be a basis 
of V . Let r < s and 

V = Wi-\ \-Wr- {Ws+1 H h Wd). 

Let c be the cone generated by Wi 's for 1 < i < d and let t be the cone 
generated by c and the vector v. 

(a) Define the semi-open cones ti by 

ti = c{wi, . . . , Wi„i) + c{wi+i, . . . ,Wd,v) for 1 < i < r, if r > 1, 
ti = c{wi, . . . , Wi^i) + c(wi+i, . . . ,Wd,v) for s + 1 < i < d, if s < d. 

Then: 

(a.l) If r > 1, t is the disjoint union of the ti for 1 < i < r. 
(a.2) t\c is the disjoint union of the ti for s + 1 < i < d. (If s = d, 
thent\c = 0.; 

(b) Forie{l,...,r}U{s + l,..., d}, let 

Ci = C{wi, . . . , Wi-i, Wi+i, . . .,Wd,v). 

Then 

r d 

M - E [^*] + E [^^] = 

i=l i=s+l 

modulo indicator functions of lower dimensional cones. 

Proof. Let us prove (a.2) first. By quotienting out by the subspace 
generated by Wi, . . . ,Ws, we can assume that s = 0, thus v = —{wi + 
. . . + Wd). In this case, we just apply Lemma 12 with Uk = Wk and 
Ud+i = V. As the cone t-d+i of Lemma 12 is just c, we obtain (a.2). 
Let us now assume that / > 1. We write 

Wi = V + Wd-\ h Ws+l - {Wr-\ h ^2)- 

We apply (a.2) with wi = v, W2 = Wd etc., and v = wi. The correspond- 
ing semi-open cones are ti = c{v, Wd, . . . , Wr+i) + ^(^,—1, . . . , W2,wi) = 
tr, ■ ■ ■ ,ir = c(f , Wd, . . . , W2) = ti, thus we obtain (a.l). 

For i G {1, . . . , r} U {s + 1, . . . , d}, the (closed) cone Cj differs from 
the semi-open cone ti by the union of some faces of lower dimension, 
therefore (b) follows from (a) and inclusion-exclusion relations. □ 
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Next, we obtain decompositions of any simplicial cone into full- 
dimensional cones which have a given edge, or a given facet. Recall 
that C{V) is the linear span of indicator functions of cones which con- 
tain straight lines. 

Proposition 14. Let wi, . . . ,Wd be a basis of V and let c be the cone 
generated by Wi 's for 1 < i < d. Let 0<r<s<d+l. 

a) Let 

V = Wi-\ \- Wr - {Ws+l -\ \-Wd). (5) 

Then we have the following relation between indicator functions of 
cones. 

r 

[c] = ^{-ly-^ [C(-Wi, . . . , -Wi_i, Wj+i, ...,Wd, v)] 

1=1 

d 

+ ^ (-1)°'"* [c(wi, . . . , Wi-i, -Wi+i, -Wd, -f )] mod C{V). 

i=s+l 

(6) 

b) Let L G V be a hyperplane. Assume that Wi E L if and only if 
r + 1 < i < s, that the Wi 's lie on one side of L for 1 < i < r and on 
the other side for s + 1 < i < d. For i G {1, . . . , r} U {s + 1, . . . , d}, 
denote by pi the projection V ^ L parallel to Wi. Then we have the 
following relation between indicator functions of cones. 

r d 

[c] = + P^{c)] - i^A-Wi) + Pi(0] mod C{V). (7) 

1=1 i=s+l 

Example 15. c = c{wi, W2), v = wi + W2- Then 

[c] = - [c{-W2, v)] + [c{wi,v)] + [C{W2, -W2, Wi)] . 

Indeed [c{w2, —W2, wi)] — [c] is the indicator function of the quadrant c{wi, —W2) 
minus that of the hne Mwi. This is also the case for [c(— -^2; ^^)] + [£(^^1, ^^)]- 
See also Figure 2. 

Example 16. We consider the 3-dimensional cone 

c = c((-l,0,0), (-1,2,0), (-1,0,3)). 

i) First let L be the subspace with basis (0, 1, 1). Formula (6) gives 

[c] = [c((-l, 0, 0), (-1, 2, 0), (0, 1, 1))] - [c((-l, 0, 0), (1, 0, -3), (0, 1, 1))] 

+ [c((-l,2,0), (-1,0,3), (0,-1,-1))]. 

ii) Now, let L be the subspace with basis (1, 0, 0), (0, 1, 1). Formula (7) gives 
[c] = [c((-l, 2, 0), (-1, 0, 0), (-5, 6, 6))] - [c((l, 0, -3), (-1, 0, 0), (-5, 6, 6))] . 
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This example will be continued in Example 26. 

Proof of Proposition 14- Let us prove a). By quotienting out by the 
subspace generated by Wr+i, ■ ■ ■ ,Ws, we can assume that r = s. We 
write (5) as 

Wr + ■ ■ ■ + Wi — V — Wd — ■ ■ ■ — Wr+l = 0. 

By applying Lemma 12 with ui = —Wr, ■ ■ ■ , Wr+i = —v,---, Ud+i = 
Wr+i, we see that the whole space V is the disjoint union of the semi- 
open cones 

= c{Wr-l,...,-V,-Wd,...,-Wr+l), 

tr = c{Wr,. . . ,102) + C{-V,-Wd,. . . ,-Wr+l), 
tr+1 = c{Wr,. . . ,Wi) + C{-Wd,. . . ,-Wr+l), 

= c{Wr,. . . ,Wi,-v) + C{-Wd-1,. . . ,-Wr+l), 

td+1 = c{Wr,...,Wi,-V,-Wd,...,-Wr+2)- 

For any basis ui, . . . , Ud, we have 

[{c{ui, ...,Uk)+ iiuk+1, . . . , Ud))] 

= (-1)'^"^ [c(mi, . . . , Mfc, -Mfc+1, . . . , -Ud)] mod C{V). 
Therefore we obtain 

[tl] = {-l)'^[c{-Wr-l, . . . ,-Wi,V,Wd, . . . ,Wr+l)] 

[tr] = (-1)'^^'' [c{Wr,. . . ,W2,V,Wd,. . . ,Wr+l)] 
[tr+l] = {-l)'^''''^[c{Wr,...,Wi,Wd,...,Wr+l)] 

[td+l] = [c{Wr,...,Wi,-V,-Wd,...,-Wr+2)]- 

Moreover these indicator functions add up to [V] = mod C{V). Thus 
we have proven a). 

Let us show that b) is just the dual of a stellar decomposition. Let 
a G V^* be such that kera = L, and such that a{wi) > for 1 < i < r 
and a{wi) < for s + 1 < i < d. Denote the dual cone by c° C V*, 

c° = eV* : {^,x) >0 for all a; G c }. 

Let (71, . . . jd) be the dual basis of (f 1, . . . , Vd). Thus c° = c(7i, . . . , 7^). 
From the stellar decomposition of the cone c° with respect to the extra 
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edge a, we see that the expression 

[c°] - [c(7i,...,7i,...,7d,a)] 

ie{l,...,r}U{s+l,...,d} 

is a hnear combination of indicator functions of lower dimensional 
cones. If a linear identity holds for indicator functions of cones, the 
same identity holds for the indicator functions of their duals, (see 
for instance [7], Corollary 2.8). Moreover, the dual of a lower di- 
mensional cone contains a line, therefore by taking indicator func- 
tions of dual cones we obtain an equahty mod CiV). (This "dual- 
ity trick" goes back to [ ]). There remains to compute the dual cone 
(c(7i, . . . , 7j, . . . , 7rf, a))°. We see easily that if {a, Wi) > 0, then 

(c(7i, . . . , 7i, . . . , 7d, a))° = R+Wi + pi{c), 

while if (a, Wi) < 0, then 

(c(7i, . . . , 7i, . . . , 7rf, a))° = M+(-u;i) + pi{c). 

Thus we have proven b). □ 

As we will see now. Proposition 14 consists of particular cases of the 
next result, Brion-Vergne decomposition. Let L be a linear subspace 
of V. Let c be a simplicial cone in V with generators Wi, . . . , Wd- We 
consider the subsets a C {1, . . . ,d} such that the vectors Wj, for j G a, 
form a supplementary basis to L in V. In other words we have a direct 
sum decomposition 

V = ^Rwj®L. 

jeer 

The corresponding projection on L is denoted by po- : V L. The set 
of these a is denoted by B{c, L). 

Remark 17. At most (^^J = 0{d^'^) such sets a exist. Thus, if 
codimL is a fixed constant ko, this is a polynomial quantity, and the 
sets a G B{c, L) can he enumerated by a straightforward algorithm in 
polynomial time. (In practice, reverse search [ ] would be used for this 
enumeration.) 

A vector ainV/L is said to be generic with respect to the pair (c, L) 
if a does not belong to any hyperplane generated by vectors among the 
projections of wi, . . . , Wd- In other words, for any a G B{c, L), we have 

a = aaj{wj mod L), 
where Oo-j- 7^ for every j G a. 
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Figure 2. Brion-Vergne decomposition of a cone c into 
cones with a face parallel to the subspace L, modulo 
cones with lines. The vectors in the quotient V j L deter- 
mine the signs eaj- 

Theorem 18 (Brion-Vergne decomposition, [' '], Theorem 1.2). Let L 
he a linear subspace of V. Let c be a full dimensional simplicial cone 
in V with generators wi, . . . ,Wd- Fix a vector a E V/L that is generic 
with respect to the pair (c, L) and belongs to the projection of c on V/L. 
For a e B{c,L), let 

a = aa^j{wj mod L). 

Let e^j be the sign of a^j and e(cr) = rijeo-^crj- Denote by C V the 
cone with edge generators 

^ajWj for i G 0", and Paiwk) for k ^ a. 

Then we have the following relation between indicator functions of 
cones. 

[c]= ^odC{V). (8) 

o-6B(c,L) 

The theorem is illustrated in Figure 2. 

Remark 19. As explained in [9], the case where dimL = 1 reduces to 
case a) of Proposition 14- Let us show that the case where codimL = 1 
reduces to case b). A vector a in V/L is generic with respect to (c, L) if 
and only if a ^ 0. Up to renumbering, we can assume that wi, . . . ,Wr 
lie on the same side of L as a, while Ws+i, ■ ■ ■ ,Wd lie on the other side. 
The set B{c, L) consists of the singletons di = {i} such that Wi ^ L, 
ie i E {1, . . . , r} U {s + 1, . . . , d} . The sign e^^^i is +1 if i E {I, ... , r} 
and —lifiE{s + l,..., d}, therefore Equation (8) becomes precisely 
Equation (7). 
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Theorem 20 (Efficient construction of a Brion-Vergne decomposi- 
tion). Fix a number k^. There exists a polynomial-time algorithm for 
the following problem. Given a rational simplicial cone c and a ratio- 
nal subspace LofV of codimension at most ko, compute a decomposi- 
tion (8). 

Proof. By Remark 17, we can enumerate the bases a G B{c, L) in 
polynomial time. 

A rational generic vector a can be constructed in polynomial time 
as follows; this technique was already used in [5]. Fix one of the bases 
a e B{c, L) as a basis of V/ L. Then the projections of Wi, . . . .Wd have 
rational coordinate vectors Wi, . . . G 0^^° with respect to this basis. 

Consider the moment curve a{t) = (l,t, t^, . . . jt^^'-^) G M^f" for t > 
0. For any a G B{c, L) and j G cr, we have by Cramer's rule that 
the coefficient a^j{t) for a = a{t) is nonzero if and only if fa,j{t) = 
det{a(t);Wi : i E a,i j) 0. The functions fa,j(t) are polynomial 
functions in t of degree at most fco(^o ~ 1); which are not identically 
zero. Hence each has at most ko{ko — 1) zeros. Thus, by computing 
the coefficients a„j(t) for M = /co(fco — 1)1-^(^,-^)1 + 1 different values 
for t, we will find one t = t that is not a zero of any /o-j. Clearly this 
search can be implemented in polynomial time. 

Then the vector a = mod L) E V/L is generic with re- 

spect to the pair (c, L), and since ai(t) > 0, it also lies in the projection 
of c on V/L. (We could also use explicit deterministic constructions of 
generic vectors of polynomial encoding size, without search, similar to 
[11, 12].) □ 

Remark 21 (Relation to variants of Barvinok's decomposition). Barvi- 
nok's method to efficiently compute the discrete generating function 
S{s -|- c)(^) of a cone in fixed dimension decomposes a simplicial cone 
into cones with smaller indices (and ultimately into unimodular cones). 
The original method by Barvinok [ ] used a primal decomposition and 
inclusion-exclusion to take care of lower-dimensional cones. The dual 
of the stellar decomposition is used in the variant of Barvinok's method 
that was presented in [ ] and later implemented in LattE [ ]; this 
avoids having to take care of lower dimensional cones. Instead of du- 
alizing, one can directly uses Formula (6), the case where dimL = 1. 
We use this primal variant of Barvinok's decomposition in our imple- 
mentation. Other primal variants were studied in [11] and [12] and 
implemented in LattE macchiato. 



3.3. Complexity of the intermediate generating function. We 
will use the following notations. 
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Definition 22. Let 

oo n 

T{t,x)^e'^-^^^-J2Bn{t)'^, (9) 

n=0 

where Bn{t) are the Bernoulli polynomials. 

Definition 23. For a and q a positive integer, we denote 

:=?Lg^J ^ {^}</ •= (nmodq) G [0,g), 

which give the unique decomposition 

[n\q + {n}q. 

Thus [aji = |_aj and {a}i = {a} will be the ordinary "floor" and 
"fractional part" of a. 

Theorem 24 (Short formula for S^{s + c)(^)). Fix a non-negative 
integer ko. There exists a polynomial time algorithm for the following 
problem. Given the following input: 

(11) a number d in unary encoding, 

(12) a simplicial cone c = c{vi, . . . ,Vd) C M*^, represented by the 
primitive vectors Vi, . . . ,Vd & l^"^ in binary encoding, 

(13) a rational subspace L C M.'^ of codimension ko, represented by 
d — ko linearly independent vectors hi, ... , ba-ko ^ binary 
encoding, 

compute the following output in binary encoding: 

(01) a finite set F, 

(02) for every 7 in F, integers a^"^^ and lattice vectors w^'^^ e Z*^ for 
i = 1, . . . ,d, 

which have the following properties. 

(1) For every 7, the family {w["'^), for i = 1, . . . ,d, is a basis o/M'^ and 
if we denote the dual basis by {rj^"'^), for i = 1, . . . ,d, then rlf'' e 1/ 
fori = l,...,ko. 

(2) For every s e M.'^, we have the following equality of meromorphic 
functions of ^: 

=e«-'^««nr({-(.w,.>},(«.-r'» ■ , (10) 

7er i=i lii=i\S)'"^j / 
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Remark 25. Consider the term corresponding to 'j G T in (10). Define 

ko 



,(7) 

i=l 

,(7) 



The point s + s^'^' is in Zw^ © ©j=fc(,+i Mtfj , formula (10) 

reads a/so 

^^(s + c,Z^)(0 = 

V «We<«'^+^'''> ^ ^ (11) 

Formula (11) zs easier to grasp, while (10) zs nsed m the program, 
where we avoid all the singular hyperplanes {^,wl'^^) = 0, for 7 G F 
and i = 1, . . . ,d, by deforming the linear form ^ into ^ + e. 

Proof. Using Theorem 18, we construct a signed decomposition of c 
into cones c"" which have a face parallel to L, so that for any s G M*^ we 
have 

[s + c] = e'^ [s + c'^] (modulo cones containing lines), (12) 



where e'^ G {±1}- Here a runs through the set B{c,L) of Theorem 18, 
which can be enumerated in polynomial time by Remark 17. This 
makes the construction of the signed decomposition a polynomial-time 
algorithm. As the valuation vanishes on cones which contain lines, 
we get 

creB{c,L) 

Each cone c'^ is given by its edge generators wf , i = 1, . . . ,d, which are 
such that (w.f ), for i = ko + 1, . . . , d) , is a. basis of L. Then we apply 
Theorem 28 [Short formula for S^^{s + c,7j'^){C,)] of [ ] with input c"" 
and / = {fco + 1, . . . , d}. See Examples 6 and 7. □ 

Example 26 (Continuation of Example 16). Again c is the cone 

c((0, -2,0), (1,-2,0), (0,-2,3)). 

Here is the output of our Maple program computing the generating functions 
S^{s + c)(^) (Formula (11)). We set s = (a, b, c) and ^ = (x, y, z). 



17 



i) Let L be the subspace with basis (0, 1, 1). In this case, the three cones 
which occur in the decomposition of Example 16 have non-unimodular pro- 
jections in V/L, so that the projections are further decomposed into uni- 
modular cones. This is why we have six terms in the expression of S^{s + c): 



S^{s + t){x,y,z) 



^ax+by+cz 



^{■ia-b+c}{ly-lz)+{a}{-x+?, z) 



;i_eif-|^)(l-e-^+3^)(y + z) 



+ 



^{a}(-x+2y)+{2a+b~c}{~ly+lz) 

;i_e-^'+2 2/)(l_e-i2'+f^)(y + z) 



{3 a—b+c}z+{—a}{x—?iz) ^—{~b+c}z+{—a}x 

~ (1 - e-^)(l - e^'-3^)(-y - z) ^ (1 - e-^)(l - e^)(-y - z) 

^{-b+c}y-{a}x ^{a}(-x+2y)-{2 a+b-c}y 

+ -71 7:^7--. 7-r, T + 



(1 - ey){l - e-^)(-?/ - z) (1 - e-^+2y)(l - e-y){-y - z) 
ii) Let L be the subspace with basis (1, 0, 0), (0, 1, 1). Then: 

S^{s + c){x,y,z) 

f p{-f'+c}(-5^+?^) 



1 , 



1 _ e-2^+yj x(-5x + 6y + 6z) 

{-b^c}{\x-z) 

+ 6^ 



1 , 



Example 27. We consider the four triangles Tj of Figure 3 which sub- 
divide the square q = [0,4] x [0,4]. Let L = M(0,1). Let us apply the 
inclusion-exclusion formula together with the valuation property of . As 
vanishes on the edges and vertices which occur in the intersections of the 
triangles, it follows that the sum 

E^=i must be equal to 5^(q). 

The example has been chosen so that the intermediate valuation 5^(q) 
can also be computed directly by elementary means. This allows us to verify 
our methods on the example. A direct computation gives 

5^(q) = + + e^" + e^- + e^"). 

Table 1 shows the output of our Maple program that implements the algo- 
rithm of the theorem for S^{Ti), z = 1, . . . , 4. The equality "^f^i S^{Ti) = 
S'^(q) is readily checked on these expressions. 
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Figure 3. The square q subdivided into 4 triangles Tj. 
Lattice subspaces parallel to the subspace L — M(0, 1) 
intersect the walls between the triangles at most in a 
point, and thus vanishes on the intersections of the 
triangles. 

Table 1 . Intermediate valuations of the square and the 
4 triangles 



p Vertices V(p) Intermediate valuation S'^(p) 



T2 

q Q - (^) ' (I) ' (o) ^ (1 + e- + e2- + e^- + e^-) 



e e _|_ 



(l-e^'+Tzy'^y (i-e^'-Tiy'^y (^i-e^+is^^j 



g3a;+§H j,3a:+4i/ ^4 a: g4 x+4 ji 

g2x+4y _ e^^^+lt^ _|_ e'^y _ e^y 



y 
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4. Ehrhart quasipolynomials for intermediate valuations 

Let us fix a rational polytope p and a weight h{x). It is well known 
("Ehrhart's Theorem") that when we dilate p by an integer t E N, 
the function t S(tp, h) is a quasi-polynomial function of t of degree 
d + M where M is the degree of /i, 

d+M 

S{tp,h)= Y,Em{p,h-t)t'^. (13) 

m=0 

The coefficients Em{p,h;t) depend only on {t}q where q is an integer 
such that qp has lattice vertices. When p is a lattice polytope, then the 
coefficients Em{p,h]t) do not depend on t, and t (-)■ S{tp,h) becomes 
a polynomial. 

E. Linke [13] has proved that the expression (13) is still valid for real 
dilations t G M. This is the generalized setting that we are going to use 
in this section. Note that if we allow real dilations, then even for lattice 
polytopes p we obtain a quasi-polynomial rather than a polynomial; the 
fractional part function {■}q with q = I appears in the expressions for 
the coefficients Em{p, h; t). For instance, the number of integers in the 
interval [0, t] = t[0, 1] is t + 1 - 

We are going to extend Ehrhart's Theorem to the intermediate val- 
uation S^{p). Indeed, one still has 

d+M 

S\tp,h) = Y,EL{p.h-t)t^ (14) 

?Ti=0 

where the coefficients E^^ip, h; t) depend only on {t}q. Note that when 
L = V, then S^{tp,h) is polynomial, and coincides with the integral 
of h on tp. 

More precisely, we will show that the Ehrhart coefficients are step- 
polynomial functions of t, in the sense of [1 1] . In the next theorem, 
we prove these results and furthermore, we show that the coefficients 
E^{p,h;t) are given by short formulas, provided that in the input 
(L,p, h), the subspace L has fixed codimension, p is a simple polytope 
given by its vertices, and the weight is a power of a linear form h{x) = 

Theorem 28. For every fixed number G N, there exists a polynomial- 
time algorithm for the following problem. 
Input: 

(Ii) a number d in unary encoding, with d > k, 
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(12) a rational subspace L C M'^ 0/ codimension ko, represented by 
d — ko linearly independent vectors hi, ... , bd-ko ^ Q'^ binary 
encoding, 

(13) a finite index set V, 

(14) a simple rational polytope p, given by the set V of its vertices, 
rational vectors Si G Q'^ in binary encoding, 

(15) a rational vector £ ^ Q"^ in binary encoding, 
(Ig) a number M G N m unary encoding. 

Output: 

(01) a finite index set T, 

(02) polynomials p'"^ G Q[ri, . . . , t^q] and numbers Q'"^ G Z, g,/'™ G 
N for 7 G r, z = 1, . . . , fco and m = 1, . . . , c? + M , 

such that fort eR we have S^{tp, {i,x)^) = where 
Emit) is given by the step-polynomial 

Emit) = E/''™({cr^k-'---'{cr%r)- 

Remark 29. Similar results hold for a more general weight h{x), as 
in [o] for integrals over a simplex. One can assume that the weight 
is given as a polynomial in a fixed number D of linear forms, h{x) = 
x), . . . , {Ed, x)), or has a fixed degree D. 

Proof of Theorem 28. As first observed in [ 1] , the path from generating 
functions to Ehrhart quasi-polynomials rehes on the following impor- 
tant property of functions 0(^) G M.(^{y*): such a function has a unique 
expansion into homogeneous rational functions 

r?i>mo 

If P is a homogeneous polynomial on V* of degree p, and D a product 
of r linear forms, then is an element in J^iiV*) homogeneous of 
degree m = p — r. For instance, |^ is homogeneous of degree 0. On this 
example, we observe that a function in Aii{V*) which has no negative 
degree terms need not be analytic. 

For a fixed t G M, the polynomial function ^ S^{tp, ^^'^j, ) is the 
term of degree M of of the holomorphic function S^{tp){C,). Let be 
the supporting cone of p at the vertex s. By Brion's theorem applied 
to the semi-rational polytope tp, we write S^(tp){^) as the sum of the 
intermediate generating functions of the cones at the vertices ts, s G V, 
of the dilated polytope tp. The crucial point is that the cone does 
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Table 2. Intermediate Ehrhart quasi-polynomials of 
the square and the 4 triangles 



p Ehrhart quasi-polynomial S^{tp, 1) 

Ti f t'-l ({4t}i)'-^{-5t}2 + i{4t}i + ^ ({-5t}2f 

T, 3t^ + (2 - 4{4t}i) t - I ({-5t}2)' + i{-5t}2 + | ({4t}if - |{4t}i 
T4 it_^ ({5t}2)2 + ^{5i}2 

q 16*2 + (4-4{4t}i)t 



not change when the polytope is dilated. 

sev 

Let gs be the smallest integer such that qgS is a lattice point. Then 
(t — is a lattice point, therefore we have 

= e*<«'^>e-«'^*>-^>S^({t},,s + c,)(0. 

Expanding e*^^'''^ in powers of t and looking for the coefficient of t\ 
we obtain 

i>o sev 

Thus we must analyze the term of ^-degree M — i of 

We use Formula (10) and consider the summands indexed by 7 G F, 
one at a time. From then on, the proof is entirely similar to the proof 
of Theorem 37 in [:'], and we omit it. □ 

Example 30 (Example 27, continued). Table 2 shows the output of our 
Maple program for the Ehrhart quasi-polynomials of the four triangles of 
Example 27 with respect to the weight h = 1. Note that the dilating pa- 
rameter t is real. We obtain 
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4 

Y^S^(tTi) = 16t'^ + At-4t{4t}i -^{-5t}2 

= 16t2 + 4t-4t{4t}i, 

using a relation between fractional parts, {x}^ — {—x}^ = {x} — {—x}, for 
the simplification. Indeed, a direct calculation gives the same answer. 

Setting aside considerations of efficient computation, we can prove 
the following theorem. 

Theorem 31 (Real Ehrhart Theorem). Let p C &e a rational poly- 
tope, h be a polynomial function of degree M , and L C'R'^ be a rational 
subspace of codimension k^. Then t h- )■ S^{tp,h) is a quasi-polynomial 
function oftE M>o- More precisely, 

d+M 

5^(ip,/i)= J]£;^(p,/i;i)i-, (16) 

m=0 

where E^{p, h; t) is a step-polynomial of the form 

Emit) = E ■ ■ ■ ' {CrO<-) for t e R>o, 

where T is a finite index set, with polynomials /^'™ G Q[?"i, • • • , t^q] 
and numbers ^7'™ ^ ^7 ?7'™ ^ for 7 G T, i = 1, . . . ,ko and m ~ 
l,...,d + M. 

Proof. The polytope p is no longer assumed to be simple, so we compute 
a triangulation of each of the vertex cones and use inclusion-exclusion 
to avoid overcounting. Using a decomposition of h into powers of linear 
forms, then the method of the previous theorem gives the result. □ 
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